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Abstract 

The Euler equation (EE) is one of the basic equations in many physical fields such as fluids, 
plasmas, condensed matter, astrophysics, oceanic and atmospheric dynamics. A symmetry group 
theorem of the (2+l)-dimensional EE is obtained via a simple direct method which is thus uti- 
lized to find exact analytical vortex and circumfluence solutions. A weak Darboux transformation 
theorem of the (2+l)-dimensional EE can be obtained for arbitrary spectral parameter from the 
general symmetry group theorem. Possible applications of the vortex and circumfluence solutions 
to tropical cyclones, especially Hurricane Katrina 2005, are demonstrated. 



PACS numbers: 47.32.-y; 47.35.-i; 92.60.-e; 02.30.1k 



I. INTRODUCTION. 



There are various important open problems in fluid physics. One of the most important 
problems is the existence and smoothness problem of the Navier-Stokes (NS) equation. The 
NS equation has been recognized as the basic equation and the very starting point of all 
problems in fluid physics [l]. Due to its importance and difficulty, it is listed as one of the 



millennium problems of the 21st century 

One of the most significant recent developments related to the above problem may be 
the discovery of Lax pairs of two- and three- dimensional Euler equations (EEs) which are 



the limit cases of the NS equation for a 



dimensional EE, the Lax pair given 



in |a, 4| i 



arge Reynolds number {3, 4]. Actually, for the two 



is weak (see Remark 2 of the next section) while 



the Lax pairs of the three-dimensional EE are strong (see Theorems 3 and 4 of which 
can be proved in a similar way as Theorem 1 of this paper). Hence, the EEs are (weak) Lax 
integrable under the meaning that they possess (weak) Lax pairs, and subsequently the NS 
equations with large Reynolds number are singular perturbations of (weak) Lax-integrable 
models. 

The (3+l)-dimensional EE 

E 1 = u} t + (u- V)u; - (u ■ V)m = 0, (1) 
u = V x u, (2) 

— * 

with V • u = is the original springboard for investigating incompressible inviscid fluid. In 
Eqs. §B) and (J2J), u = {tu 1 , uj 2 , ^3} is the vorticity and u = {ui, u 2 , %} is the velocity of 
the fluid. 

In (2+l)-dimensional case, the EE has the form of 

E = U) t + [iP, LO\ = 0, LO = lfj XX + ifjyy, (3) 

where the velocity u = {u\, u 2 } is determined by the stream function if) through 

Ul = -Ipy, U 2 = lfj x (4) 

and the Jacobian operator (or, namely, commutator) [A, B] is defined as 

[A, B] = A x B y - B x A y . (5) 



It is known that the EEs are important not only in fluid physics [6j but also in many 



other physical fields such as plasma physics [7|, oceanography 



161 



, atmospheric dynamics 



9j, superfluid and superconductivity Il0|], cosmography and astrophysics fll|. statistical 
physics field and particle physics]l3[] and condensed matter including Bose-Einstein 
condensation [y|], crystal liquid 15] and liquid metallic hydrogen 16], etc.. 

As a beginning point of various physical problems, the EEs have been studied extensively 
and intensively, which is manifested by a large number of related papers on EEs in the 
literature. For instance, a lot of exact analytical solutions of the EEs have been presented, 
some of which can be found in the classical book of H. Lamb In the authors studied 
the planar rotational flows of an ideal fluid and the addressing method was developed to 



obtain exact solutions of the EEs in 



191 ] . In addition, some types of exact solutions were 



obtained via a Backlund transformation in 20[ . However, rather few exact analytic solutions 
of the EEs have been obtained from the (weak) Lax pair since it was revealed by Charles Li 
{3I around five years ago. A special type of Darboux transformation (DT) with zero spectral 
parameter for the (2+l)-dimensional EE was shown in and some types of DTs (or weak 
DTs) with nonzero spectral parameter (s) for both (2+1)- and (3+l)-dimensional EEs were 
presented in our unpublished paper 5]. 

Lie group theory is one of the most effective methods of seeking exact and analytic 
solutions of physical systems. However, even for a mathematician, it is still rather difficult to 
find a symmetry group, especially non-Lie and non-local symmetry groups. So for physicists, 
it would be more significant and meaningful to establish a simple method to obtain more 
general symmetry groups of nonlinear systems without using complicated group theory. 

To our knowledge, there is little exact analytic understanding of the vortices and circum- 
fluence, although they are most general observations in some physical fields; in particular, 
very rich vortex structures exist in fluid systems. In fact, if one could find the full symmetry 
groups of the EEs, then many kinds of exact vortex and circumfluence solutions could be 
generated from some simple trivial solutions. 

This paper is an enlarged version of our earlier, unpublished paper 5|]. In section II, 
we first establish a simple direct method to find a general group transformation theorem 
for the (2+l)-dimensional EE, then utilize the theorem in some special cases to obtain 
some solution theorems which lead to a quite general symmetric vortex solution with some 
arbitrary functions. The applications of the exact vortices and circumfluence solutions are 



given in section III. It is indicated that the solutions can explain the tropical cyclone (TC) 
eye, the track, and the relation between the track and the background wind. The TC tracks 
can thus be predicted by the relation. In section IV, beginning with a general symmetry 
group theorem, the DT in with zero spectral parameter is extended to that with arbitrary 
spectral parameter. The last section is a short summary and discussion. 



II. SPACE-TIME TRANSFORMATION GROUP OF THE TWO-DIMENSIONAL 
EE. 

In the traditional theory, to find the Lie symmetry group of a given nonlinear physical 
system, one has to first find its Lie symmetry algebra and then use the Lie's first fundamental 
theorem to solve an "initial" problem. If one utilizes the standard Lie group theory to study 
the symmetry group of the two-dimensional EE, it is easy to find that the only possible 
symmetry transformations are the arbitrary time-dependent space and stream translations, 
constant time translation, space rotation and scaling j^ll ]. 

Recently, for simplicity and finding more general symmetry groups, some types of new 
simple direct method without the use of a ny g roup theory have been established for both 



Lax-integrable 22J and non-Lax-integrable 23J models. 



For the two-dimensional EE ([3] 
Theor 
Lax pair 



, we have the following (weak) Lax pair theorem. 
Theorem 1 (Lax pair theorem [S]). The (2+l)-dimensional EE ([3]) possesses the weak 



u x (j)y - Uy(j) x = \<f>, (6) 

4>t + 1px(t)y ~ 4>y(f>x = 0, (7) 

with the spectral parameter A. 

Proof. To prove the theorem, we rewrite ([6]) and (J7J) as 

L(p = 0, L = [u, ■] - A, (8) 

M0 = O, M = d t + [^, •]. (9) 

It is straightforward that the compatibility condition of Eqs. (JED and (El), LM — ML = 0, 

reads 

LM — ML = -[u t , ■} - [V, k ■]] + k fofr, ■]] = 0. (10) 



Using the Jacobian identity for the commutator [•, •] defined by (J5]) 



[A, [B, C)) + [B, [C, A]} + [C, [A, B]] = 0, 

Eq. ( TTOT) becomes 

[ut+lij, u>], -] = 0. (11) 

The theorem is proven. □ 



Remark 1. The theorem was proved in a slightly weak way in [4], where the compatibility 
condition of Eqs. and (J7J) was 

[u t + [iji, u), 0] = O, (12) 

with the requirement that was just the spectral function. However, in our new proof 
procedure, the spectral function in Eq. (1121) can be replaced by any arbitrary function. 
Remark 2. In Theorem 1, the Lax pair is termed weak because starting from the Lax pair, 
we can only prove Eq. ( flT|) instead of the EE ([3]) itself. For instance, Eq. ( fTTj) is true for 

oj t + [ip, u)] = c(t) 

with c(t) being an arbitrary function of t. Therefore, all the conclusions obtained from the 
Lax pair have to be treated carefully by substituting the final results to the original EE to 
rule out the additional freedoms. 

From Theorem 1, we know that the (2+l)-dimensional Euler equation is weak Lax inte- 



grable. So we can apply the new direct method developed in [22] to find some complicated 
exact solutions from some simple special trivial ones after ruling out the ambiguity men- 
tioned in remark 2. 



Using the method in [221 ] . we have the following transformation theorem: 

Theorem 2. (Group Theorem). If {uj'(x, y, t), ifi'(x, y, t), (f>'(x, y, t)} is a known solution 

of the two-dimensional EE ([3]) and its Lax pair ([6]) and ([7j) with the spectral parameter A', 

{u, ip, 0} with 



= exp(#)0'(£, 77, r) = exp(#)0' (13) 
is a solution of Eq. fTTT]) and its Lax pair with the spectral parameter A, if and only if the 



following three conditions are satisfied: 

{([t,u]^ + [u,rj\)\'- Ao4+ [uJ,g]oj' t }(t>' 

+ { [r, u]u/ T + [w, vW v + k e]^} ^ = °> ( 14 ) 
{ fa + 77] - (r t + [V, r])^] A' + ( 5t + [V>, <?]K} 0' 

+ {(r t + hMK + ( Vt + tyMWr, + (6 + [MW(] 4 = °> ( 15 ) 

^ = Vw + Vw> ( 16 ) 

where the arguments {x, y, t} of the functions u', ip' and <f>' have been transformed to 
{£> T }> an< i £> Vi T an d 5 are functions of {x, y, t}. 

Proof. Because {u)'(x, y, t), ip'(x, y, t), <f/(x, y, £)} is a solution of the EE and its Lax 
pair with the spectral parameter A', then {u/(£, r], r), Vi T )> 0'(£> ?7> r )} satisfies 

^;- w ;4 = a>', (17) 

0; + ^;-^ = o, (is) 

and 

w ; + - ^ = 0. (19) 

Substituting Eq. ( fT3l) into Eqs. © and (EJ), we have 



[u, + [w, ^ + k r]# + ([w, </] - A)<// = 0, (20) 

fa + [</>, e])^ + fa + + fa + k r])# + fa + [V, = 0. (21) 

Applying Eqs. (jTTj) . ffTBl and (H2J to Eqs. (J2DD and (ETJ) by ruling out the quantities <// and 
<f)' v yields Eqs. (Ill]) and (fl5l). 



It is noted that Eq. (1161) in Theorem 2 is only the definition equation of the vorticity. 
Theorem 2 is proven. □ 

From Theorem 2, we have only three determinant equations for six undetermined func- 
tions £, T], t, ip, uj and g, which means that the determinant equation system (|T6|) is 
underdetermined. Therefore, there exist abundant interesting exact solutions. Here we 
consider two special interesting cases of Theorem 2. 

Corollary 1. If if)'(x, y, t) is a solution of the Poisson equation 

u = < x + ij' yy (22) 



with a constant Uq, then {u, if)} is a solution of (fTTj) if the following three conditions hold: 

[T,w]Ve + M = o, (23) 

7fc + hM]-(r t + hM)^ = 0, (24) 
w = ^ + (25) 

where ■0' = i(/(x, y, t) has been redefined as tp'(^, rj, r). 

Proof. It is clear that the EE ([3]) [and then Eq. (fl~T|) ] possesses a trivial constant vorticity 
solution {a;', tp'} = {uj , tp'} with if;' being a solution of the Poisson equation. Substituting 
to' = ojq = const, into Theorem 2 results in the Corollary 1 at once. □ 

Corollary 2. If {u>'(x, y, t), if)'(x, y, t)} is a known solution of the two-dimensional EE 
([3]), then {u, -0} with the conditions 

[r, u]<J T + [u, 7]W V + [w, £]<4 = 0, (26) 

(r t + r]K + + [V, + (6 + £])<4 = °> ( 27 ) 

W = ^XX + Ipyy (28) 

is a solution of Eq. (ITTj) . where the arguments {x, ?/, t} of the functions a;' and ^// have 
been transformed to {£, 77, r}, and £, 77 and r are functions of {x, y, £}. 

Corollary 2 can be readily obtained from Theorem 2 by taking A' = A = g = 0. 

Remark 3. Corollary 1 and Corollary 2 are independent of the Lax pair though they are 
derived by means of the Lax pair. 

By solving Corollary 2, we can get the following theorem. 

Theorem 3 (Solution theorem). The (2+l)-dimensional EE possesses a special solution 
{u>, ip} with 

u = F(f(x, y, t)) = F, (29) 

V = G(f(x, y, t),t)- [ X f *(*> r > *} dz = G + h, (30) 

J fr{z, r, t) 

where / = f(x, y, t), h = h(x, y, t), F(f) and G(f, t) are functions of the indicated 
variables, the variable r = r(x, y, z, t) is determined by 

f(z, r, t) = f(x, y, t) 

and the functions F, G and / (h) are linked by the following constrained condition 

F = G f f{fl + fy) + <?,(/„ + fyy) + h XX + hyy. (31) 



Proof. After rewriting Eqs. (126|) and ([27j) as 

[w, a/] = 0, (26') 

J t + a;'] = 0, (27') 
it is not difficult to find that the general solution of Eq. (1261) [i.e. Eq. (26')] is 

oo = F(uj', t). (32) 

Though u'(x, y, t) should be an exact known solution of the EE, u/(£, rj, r) = / can still 
be considered as an arbitrary function of {x, y, t} due to the fact that £, 77 and r are all 
undetermined arbitrary functions of {x, y, t}. Then Eq. (132]) becomes 

w = F(f, t). (33) 

The general solution of Eq. (27') [or Eq. ([27]) ] is rightly Eq. (130]) . while Eq. ([31]) is just the 
direct substitution of Eqs. ([33]) and flU to Eq. (1281). 

Finally, to rule out the ambiguity brought by the weak Lax pair by substituting Eqs. 
([33]) and ([30]) into Eq. ([3]), one can find that Eq. ( ]33j) with Eq. ([30]) is really a solution of 
the EE ([3]) only if F(f,t) = F(f). Theorem 3 is proven. □ 

Because of the arbitrary function /, we can obtain many physically interesting solutions 
from Theorem 3. For instance, if the arbitrary function / is assumed to be the form 

f = (x- x ) 2 + (y- y ) 2 + h = r + h , (34) 

where x , y and h are all arbitrary functions of t, then we readily have the following special 
solution theorem. 

Theorem 4 (Special solution theorem). The (2+l)-dimensional EE ([3]) possesses an exact 
solution 

i> = Votx - x ot y + Ft In r + F 2 - -h ot tan h - / dr, (35) 

2 y- yo 4 J r 

uo = F r (r + h ), (36) 

where xo, yo, h , F± and F 2 are arbitrary functions of t, and F = F{r + h ) is an arbitrary 
function of r + ho- 

The intrusion of many arbitrary functions into the exact solution ( [35]) allows us to find 
various vortex and circumfluence structures by selecting them in different ways. 
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FIG. 1: The structure of the singular vortex expressed by (|37p with the parameters (|38p . The 
length of the arrow stands for the strength of the velocity field and the values from inside to 
outside are 16/3, 8/3, 16/9, 4/3, 16/15, 8/9, 16/21 and 2/3, respectively. 

In the solution (!35j) . the first two terms 

y ot x - x ot y 

represent the background wind (induced flow) with the time-dependent velocity field 

u = {x ot , y ot }. 

The third term ( independent ) 

Filnr, (37) 
corresponds to a time-dependent singular vortex. The detailed velocity field with 



Fi = l, 



x = y = 



(38) 



is shown in Fig. 1. All the quantities used in the figures of this paper are dimensionless 
except for the special indication in Fig. 8. 

The fourth term F2 is trivial because of the existence of the time-dependent translation 
freedom when one introduces the potential of the velocity — i.e., the stream function. 

The fifth term (/i rdependent) 



1 _! X - X 

-h t tan 

2 y-yo 



(39) 



FIG. 2: (a). The structure of the hole expressed by ([39]) with xq = yo = 0, ho = t. (b). The source 
structure of (|39p with xq = yo = and ho = —t. The length of the arrows expresses the strength 
of the velocity field and the values from inside to outside are 0.71, 0.32, 0.20, 0.14, 0.11, 0.09, and 
0.08, respectively, both for (a) and (b). 

is related to a hole [Fig. 2(a)] or a source [Fig. 2(b)]. 
The last term of Eq. (1351) 

1 fF(r + h ) 



dr, (40) 

is the most interesting because it is related to abundant vortex structures due to the arbi- 
trariness of the function F. Here are some special examples based on the different selections 
of the arbitrary function. 

(i) Lump-type vortices. If the function F{r) is a rational solution of r, 



ZloW* _Pi 



f(rH S^ s wy (41) 

with the conditions ^ and Q(r) ^ for all r > 0, then the solution (j4TJl) becomes an 
analytical lump- type vortex and/or circumfluence solution for the velocity field. Figure 3 
displays a special lump-type vortex structure of the velocity field described by ( HOI) with 

F(r) = -^-j, x = y = 0. (42) 

(ii) Dromion-type vortices. When the function F(r) is fixed as a rational function of r 



1.5 
1 

0.5 

y oj 

-0.5 
-1 
-1.54 



-1.5 



/ / ///-- -v\\\ 

I //// \\\\ 



I 



li 

WW- 

\\\\\- 
\ \\\-»- 



\ \ 
II 
"ft, 
-•"/// 

— W// 



\ \ \ 

' 1 I 



/ / I 
/ / / 



-0.5 



0.5 



1.5 



FIG. 3: A typical lump-type vortex expressed by (|40p with (|42p . The strength of the velocity field 
is expressed by the length of the arrows and the values from inside to outside are 0.44, 0.87, 1.26, 
1.52, 1.60, 1.48, 1.26, 1.01, 0.80, 0.62, 0.49, 0.39, 0.31, 0.25, 0.21, 0.17, and 0.14, respectively. 

multiplied by an exponentially decaying factor — for instance, 



(43) 

with arbitrary constants c^, hi and c — then PUI) turns into an analytical dromion-type vortex 
and/or circumfluence solution. Figure 4 exhibits a particular dromion-type vortex structure 
of Eq. (40) with 



F(r) = rexp(— r), x = y = 0. 



(44) 



iii) Rin g so litons and circumfluence. Recently, some kinds of ring soliton solutions were 



discovered 



24 



251 ]. It is interesting that the basin and plateau types of ring solitons may 



be responsible for the circumfluence solution for fluid systems described by the EE. For 
instance, if F(r) is assumed to have the property 

d l F(r) 

— tV =0, i = 0, 1, n, 

for n > 2, then (1401) expresses the circumfluence for the velocity field and the basin- or 
plateau-type ring soliton for the stream function. Figure 5(a) exhibits a special picture with 



F(r) 



-4rV 



(45) 
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FIG. 4: A typical dromion-type vortex expressed by (|40p with Eqs. (|44p . The strength of the 
velocity field is expressed by the length of the arrows and the values from inside to outside are 
0.06, 0.12, 0.17, 0.20, 0.21, 0.20, 0.19, 0.17, 0.14, 0.11, 0.08, 0.06, 0.04, 0.032, and 0.02, respectively. 

of the circumfluence structure for the velocity field, Fig. 5(b) displays the corresponding 
basin-type ring soliton shape for the stream function ip, and Fig. 5(c) shows the structure 
of the vorticity. 



III. APPLICATIONS TO HURRICANE KATRINA 2005. 

It is demonstrated in the last section that the exact solutions fl35l)- (l36l) have quite rich 
structures. Due to the richness of the solution structures and wide applications of the vortex 
in various fields such as fluids, plasma, oceanic and atmospheric dynamics, cosmography, 
assies, co— -U., etc. fl-Q , our results may be applied in all these fields. 
For instance, in oceanic and atmospheric dynamics, the analytical solution ( 1351) can be used 
to approximately describe TCs which possess increasing destructiveness over the past 30 
years 26j. The relatively tranquil part, the center of the circumfluence shown in Fig. 5(a) 



is responsible for the TC eye 



27] 



To describe different types of vortexes, one may select different types of function F(r). 
To qualitatively and even quantitatively characterize TCs, we may require that F(r) have 
the form 

F(r) = ±a 2 r 1+fe V c2 ^, (46) 
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FIG. 5: (a) A field plot of the circumfluence (|40p with Eq. (|45p for the velocity {u, v}. The 
strength of the velocity field is expressed by the length of the arrows, and the values from inside 
to outside are 0.001, 0.027, 0.131, 0.295, 0.400, 0.370, 0.248, 0.125, 0.049, and 0.015, respectively, 
(b) The corresponding basin- type ring soliton for the stream function tjj related to (a), (c) The 
corresponding bow-type ring soliton for the vorticity ui. 

with constants a, b and c. In Eq. ( )46l) . the signs "+" and "— " dictate the TCs of the 
northern and southern hemisphere, respectively. The constants a, b and c are responsible 
for the strength, the size of the TCs eye, and the width of the TC. 



FIG. 6: (Color online) The satellite image of Hurricane Katrina 2005 at 14:15, August 29, 2005, 
Coordinated Universal Time. 



The corresponding stream function related to the selection (06]) reads 

1 



V = Votx - x ot y + fjlt^ r-^e-^ { (2b 2 + 1) 



A b2 b 2 T(b 2 )T [b 2 + ^ 



-V^T(2b 2 + 1, c 2 v^)] e c2 ^ - r b2+1 2 C 4b2+2 } , (47) 

where T(z) and T(a, z) are the usual Gamma and incomplete Gamma functions, respectively. 
For more concreteness, we take Hurricane Katrina 2005 as an illustration. 



Figure 6 is the satellite image downloaded from the web |http: //www.katrina.noaa.gov/ 



satellite/satellite. html [28( for Hurricane Katrina 2005 at 14:15, August 29, 2005, Coordi- 
nated Universal Time (UTC). 

To fix the constants a, b and c in Eq. (146p for Hurricane Katrina 2005 shown in Fig. 6, 
we need know the strength (the maximum wind speed), the eye size, and the width of TC 
Katrina at 14:15, August 29, 2005 UTC. The strength of Katrina can be found in several 



websites. The data of Table I are downloaded from 



29|. 



Table I. Data of Hurricane Katrina downloaded from 



29] 
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From Table I, we know that the maximum wind speed of the Katrina 2005 at 14:15, 
August 29, 2005 UTC is about 130 mph (miles per hour) — i.e., 



^max ~ 130 mph ~ 200 km ph ~ 2 degree ph, v = yip' 2 , + ip 2 - (48) 

Comparing Katrina's satellite image shown in Fig. 6 with the map of New Orleans — say 
the map shown in Fig. 7 downloaded from 30J — one can estimate that the eye size (E) is 
about 

E <~ 1 degree ~ 100 km 
and the width (W) of the hurricane is approximately 

W fa 10 degree fa 1000 km 

for Katrina at 14:15, August 29, 2005 UTC. Using these data, we can find that the stream 
function of Katrina 2005 near the time at 14:15, August 29, 2005 UTC can be approximately 
described by 

^Katrina ~ Vot x ~ x otV ~ 4 ( 2 + 2 Vr + r) exp(-Vr), (r = (x - x ) 2 + (y - y ) 2 ), (49) 
which corresponds to the parameter selections 

a 2 ~ 8, b 2 ~ 0.5, c 2 ~ 1, 



in Eq. (j4bj) . In the real case, the quantities a, b and c should be time dependent. So the 
description here is only approximate because it is only a solution of the EE instead of the 
NS equation. 

If the strength v max , the size of the hurricane eye E, and the width W are assumed to 
have some errors, 

w max fa 130 ± 10 mph, E fa 1 ± 0.2 degree, W fa 10 ± 2 degree, (50) 

then the parameters a, b and c in Eq. (1461) or (|47p have the ranges 

a 2 fa 6.75 ~ 9.35, b 2 fa 0.25 ~ 1.25, c 2 fa 0.9 ~ 1.2. (51) 

In (I49p . xq = xo(t) and = ?/o(^) can be obtained from the data in Table I. According 
to Table I, we can find the theoretical fit of hurricane Katrina 2005 from 22:00, August 25, 
2005 to 15:00, August 30, 2005 can be approximately described by 

x = 0.00022t 2 - 0.14* - 79, y = 0.00073t 2 - 0.07t + 26 (52) 



Hurricane Katrina 



21:00 Tub August 23, 2009 to 21:00 Wed August 31, 2005 UTC 




FIG. 7: (Color online) The map with the longitude and latitude degree coordinates near the New 
Orleans and the real track of Hurricane Katrina 2005 from 21:00 Tues., August 23, 2005, to 21:00 
Wednes., August 31, 2005 UTC. 

before 21:00, August 27, 2005 and 

x = 0.0029* 2 - 0.5* - 68, y = 0.0023* 2 - 0.2* + 29 (53) 

after 21:00, August 27, 2005. In Eqs. (I52j) and (153}) . the units of x , y and t are longitude 
degree, latitude degree, and hour, respectively, while the initial time t = is taken as 22:00, 
August 25, 2005 UTC. 

Remark 4- If we fit the track only for the [Longitude, Latitude] positions, we may get 
a better fit without using any switch point. However, if we fit the track not only for the 
positions but also for times, we have to select some switch points. Physically speaking, when 
we use a parabolic line such as (152]) to fit the track of a TC, we have to assume that the 
TC moves under a constant force during the fit time period. The necessary selections of 
the switch points are caused by the fact that the driven force of the TC is time dependent. 
Here we find that if we select 21:00, August 27, 2005 as a switch point, then Eqs. (15^1) and 
(|53|) can fit the track quit well (the square error [see later, Eq. (156]) ] becomes smallest). 
This means that the TC is approximately driven by two constant forces before and after 
the turning time, respectively. Actually, approximately speaking, after this switch point, the 
TC becomes stronger and stronger (see Fig. 7 and/or Table I). 

Figure 8 describes the velocity field of Katrina at 14:15, August 29, 2005 UTC (t = 88.15) 
when the stream function is given by Eqs. ( ]49i) with (|53l) . 

In addition, the solution (|35]) - (|36]) also provides a relation between the TC track given 
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Longitude 

FIG. 8: The field and density plot for the velocity field Hurricane Katrina 2005 at 14:15/29/08 
described by Eqs. (|49p with f)53[) . The strength of the velocity field is expressed by the length of 
the arrows and the values from left to right at y = are 0.14, 0.30, 0.53, 0.82, 1.17, 1.55, 1.86, 
1.95, 1.60, 0.69, 0.18, 0.89, 1.87, 2.33, 2.31, 2.03, 1.66, 1.30, 0.99, 0.75, and 0.58 (degree per hour), 
respectively. 

by {xq, yo} and the strength of the background wind (steering flow). 

The stream function of the steering flow, ip s , can be obtained by eliminating the TC term 
(vortex term) in Eq. (1351) with F\ = F<i = ho = by setting F = and then the velocity 
field flow u s of the background wind reads 

~u s = {-ipsy, i>sx] = {x ot , Vat}- (54) 

This fact implies that once the background wind, or the large-scale steering flow in the 
upper air, is known then the motion of the hurricane center can be obtained. 

Inversely, if the motion of the hurricane center is known then the steering flow will be 
obtained at the same time. Therefore, if the position {xq, yo} of the hurricane center is 
determined, in a not very long time (say, shorter than one day), one can consider that the 
velocity of the TC will approximately keep the latest known velocity and then the TC's new 
position {xi(ti), t/i(ti)} at time t x can be determined by using 

xi{h) = x (t ) + x ot {t ){ti - t ), y x (t x ) = y (t ) + yot{to){ti - t ). (55) 

The concrete steps to predict the track and position of a TC are as follows. 



(i) Get the original known position data of a TC from professional meteorologic web site. 
The concrete position data of a happening TC, say, Katrina, can be read off from some web 
sites, say, 



29 



3l| . which are given by some international satellites and updated every six 
hours and usually three hours (or every hour) close to the landing time. 

(ii) Take the coordinate of the fit track. From the web site we can get the position 
described by longitude and latitude. Because the TCs happen in a quite small area compared 
to the whole Earth, the fit curve can be taken in a two-dimensional plane. To simplify, the 
longitude and latitude are defined as X axis Y axes, respectively. The first time recorded 
on the web site is set as initial time, and the following times are added in order by the time 
interval. 

(iii) Fit the function curve and forecast the track and position. From the first few known 
positions, it is easy to calculate the fit curve which is the function of time t. Usually we 
can take it possesses the polynomial forms of the time t, say, {X = x (t) = J2iLo a ^ N ^ = 
yo(t) = YliLo M 1 }- in this paper, we take N = 2. Finally we should minimize the square 
error, A, among the fit track and the real track 

A =Y1 l>o(*i) - x if + (l/o(*i) - Vjf] (56) 
j=m 

by fixing the constants and bi, where {xj, yj} and {xo(tj), yo(tj)} are the real and fit 
positions of the hurricane center at time tj, ni, ri\ + 1, n 2 are related to the points 
used to fit the theoretical track {x (t), yo(t)}, t = t n2 corresponds to the time to make the 
further prediction. Usually, we take t n , 2 — t ni ~ 24 (hours) that means the earlier history 
can be neglected to the hurricane track. 

Based on the above descriptions, we can use first few known position data of the hurricane 
center to predict the possible position of the hurricane some hours later. 

Figure 9 displays an example on the TC track [Fig. 9(a)] and the related background 
wind field [Figs. 9(b)-9(f)]. The zigzag line in Fig. 9(a) is the real track of Katrina 2005 
from 21:00/25/08 to 15:00/30/08 (data are read off from 2j|), and the solid line is our fit 
(by using all the data in Table I) given by Eqs. (1521 and (|53|) . Figures. 9(b)-9(f) reveal the 
corresponding steering flows at five different times. It is shown that the background wind 
leads to the change of the direction of the TC track. According to the relation between the 
track and the steering flow, we can predict the TC track by using several beginning data. 
The cross points in Fig. 9(a) are our predicted track 6 h before the real one. The same idea 
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FIG. 9: (a). The TA'gL&'g line stands for the real TC track of Hurricane Katrina 2005 from 
21:00/25/08 to 15:00/30/08 for time and from [79.9,26.2] to [87.5,36.3] for [Longitude, Latitude] 
position, the smooth line is the fit track by using all the data in Table I and the cross points 
express our predicted track six hours before the real time, (b)-(f). The corresponding steering 
flows at 23:00/25/08, 21:00/27/08, 15:00/28/08, 05:00/29/08 and 00:00/30/08 with the strengths 
1.33, 1.92, 2.92, 2.67, and 1.08 (degree per hour), respectively. 



has been applied to typhoon Chanchu 2006 33j and Hurricane Andrew 1992 



a. 



IV. DARBOUX TRANSFORMATION OF THE (2+l)-DIMENSIONAL EE 

In Sec. II, we have established a general group theorem (Theorem 2) for the (2+1)- 
dimensional EE, which can yield various solutions. Furthermore, two solution corollaries on 
the (2+l)-dimensional EE are obtained by utilizing particular seed solutions. 



It is noted that the weak DT theorems given in 



34| are special cases of the general 



group Theorem 2. In [3(, Li found a (weak) DT of the EE ([3]) with the Lax pair ([6]) and 
(J7J) for a zero spectral parameter. In 5|, the weak DT was extended to a general nonzero 
spectral parameter and many kinds of exact solutions including the solitary, Rossby, conoid, 
and Bessel waves were obtained subsquently. 

Here we derive the weak DT theorem directly from the general group Theorem 2. 

Theorem 5 (Weak DT theorem). If {a/, if)', <f>'} is a solution of the (2+l)-dimensional 
EE 03]) and its Lax pair ([6]) and ([7]) with the spectral parameter A', g(f) being an arbitrary 



function of / which is a given spectral function of (jSJ) and ([7]) under the spectral parameter 
A , then 

{u, if), </>} = W + q, ifj'+p, exp(g)cf>'} (57) 

with the spectral parameter A is a solution of the weak Lax pair (jBJ) and (J7j) and then Eq. 
( Hip where p and q are determined by 

q = Pxx + Pyy (58) 

[p, M=0, ( 59 ) 
[g, Ml+A'-A + A o /<// = 0. (60) 

Proof. Known from the proof of the group Theorem 2, Eqs. (114j) and (115j) are equivalent 
to Eqs. (EQ]) and (EE]). Taking f = x, ry = y and r = i in Eqs. (EDD and (J2U), we have 

[w, + (/]-A)0' = O, (20') 

[V>, <{>'} + <t>t + (9t + &, g\W = 0- (210 

Substituting Eq. (1ST)) into Eq. (20'), we have 

[u>, cf>'] + [q, <f>'] + (W, #] + [<?> d-A)^ 
= A'0' + [g, 0'] + ^A o /0' + ([<?, s] - A)0' 
= {A / -A + Ao/(? / + [?, ln0']}0' = O. 

Equation ( 160)) is proven. 

Substituting Eq. ([57} into Eq. (21') yields 

[if)\ <P'} + [p, 0'] + 0; + (g t + [if)\ «?] + [p, <?])0' 

= [p, 01 + [p, 5#' = [P, ln^' + = [p, ln0]0' = 0. 

Equation (1591) is proven, and Eq. fl58l) is a direct result from the definition equation of the 
vorticity. Theorem 5 is proven. □ 

It is interesting that if all the parameters A, Ao and A' are zero and the vorticity of the 
seed solution is not a constant, the Backlund transformation 

uj = oj' + q, if) = if)' + p 



with Eqs. (1591) and (!60|) is equivalent to what were obtained by Li [3] . To see it more clearly, 
one can write Eqs. fl59l) and fl60|) in the alternative forms by eliminating <fi' y via the Lax pair 
© and (JZD, 

AH - [A + Ao/^/](w' + q) x + W, q)[\nct>} x = 0, (61) 
[u' + q, p]\hy(j>) x + \ , p t = 0. (62) 

The equivalent forms of Eqs. ( l61i) and ( 1621) can also be obtained directly from Eqs. ( |T4l) 
and (|T5|) by setting £ = x, 77 = t/ and r = t. 

Remark 5. If the seed solution has a constant vorticity, the equation systems ( |59|) - (l62|) 
are completely not equivalent. Actually, when one takes a constant vorticity as a seed for 
the zero spectral parameters, nothing can be obtained from Eqs. (16TT) and (1621) . However, 
one can really find some nontrivial solutions from Eq. ( IBTfl) with a constant vorticity seed. 



In 34| , the weak DT theorem has been used to obtain some types of exact solutions such as 
the solitary waves, the conoid periodic waves, the Rossby waves, and many kinds of Bessel 
waves. Here we will not discuss them further. 



V. SUMMARY AND DISCUSSION. 

The analytical and exact forms of the vortices and circumfluence of the two-dimensional 
fluid are studied by means of the general symmetry group Theorem 2 of the (2+1)- 
dimensional EE. Some solution theorems for the (2+l)-dimensional EE are obtained from the 
group theorem by taking special seed solutions. A special weak DT of the (2+l)-dimensional 
EE is also obtained from the general group theorem. 

The special solution Theorem 4 gives a quite general exact explicit solution which covers 
many kinds of possible vortices and circumfluence such as the lump-type vortices, dromion- 
type vortices, ring solitons, etc. The vortex and circumfluence solutions may have applica- 



tions in various physical fields mentioned in the Introduction and Refs. p]-[la|. Particularly, 
they can qualitatively explain some fundamental problems of TCs such as their eye, track, 
and the relation between the track and the background wind, and the relation can be used 
to predict well the TC tracks. Hurricanes and/or typhoons have tremendously and increas- 
ingly caused destruction of our world. The method in this paper provides a possible way to 
understand and study similar disasters intensively. As an original study in this aspect, some 



introductory analysis by means of our method of Hurricane Katrina 2005 (which almost 
completely destroyed a whole city, New Orleans) are presented. The technological observa- 
tions a nd p henomenological discussions on Hurricane Katrina 2005 can be found in many 
papers 32|. In this paper, an approximate analytical expression for the (2+1) -dimensional 
stream function of Katrina 2005 is obtained. The expression is an exact solution of the 
(2+l)-dimensional EE and includes some messages including the eye size, the hurricane 
size, the strength, the relation between the hurricane center and the steering flow, etc. The 
relation is also used to predict the track of the hurricane. 

The discovery of the general group theorem may lead to the discovery of various inter- 
esting exact solutions which can be applied to many real physical fields. This paper is just 
the beginning study in this aspect. There are various important problems should be studied 
further. For instance, the possible solutions from the general group Theorem 2 are only 
discussed in three very special cases: (i) the constant vorticity seed (Corollary 1), (ii) the 
zero spectral parameters without the gauge transformation (Corollary 2, Theorems 3 and 
4), and (iii) the pure weak DT case (Theorem 5). 

In this paper, we only discuss the (2+l)-dimensional EE. Two types of Lax pairs of the 
(3+l)-dimensional EE have also been given in [4j and some special DTs of these Lax pairs 
have also been given in [5j. However, these DTs have not yet been utilized to find exact 
solutions of the (3+l)-dimensional EE. Furthermore, the corresponding symmetry groups 
similar to that of the (2+l)-dimensional EE given in this paper have not yet been discussed. 

The Lax pair and then the DT found in this paper have only weak meaning. Whether 
the (2+l)-dimensional EE is integrable under some stronger meanings [similar to those of 
(3+l)-dimensional EE] is still open. 

The more general applications of the vortex solutions given in this paper both in atmo- 
spheric dynamics and in other physical fields deserve more investigations. Especially, to 
describe the hurricane more effectively and accurately, some other important factors such 
as the Coriolis force and the viscosity of the fluid must be considered. 

Because of the importance of the EEs and the NS system and their wide applications, 
the models and all the problems mentioned above are worthy of further study. 
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